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ABSTRACT 

We use the lag-luminosity relation to calculate self-consistently the redshifts, 
apparent peak bolometric luminosities Lb, and isotropic energies for a large 
sample of BATSE bursts. We consider two different forms of the lag-luminosity 
relation; for both forms the median redshift for our burst database is 1.6. We 
model the resulting sample of burst energies with power law and Gaussian dis- 
tributions, both of which are reasonable models. The power law model has an 
index of a = 1.76 ± 0.05 (95% confidence) as opposed to the index of a = 2 
predicted by the simple universal jet profile model; however, reasonable refine- 
ments to this model permit much greater flexibility in reconciling predicted and 
observed energy distributions. 

Subject headings: gamma-rays: bursts 


1. Introduction 

The major breakthroughs in the study of gamma-ray bursts of the past six years — 
most if not all bursts are cosmological, bursts are not standard candles, the fireballs are 
beamed, many bursts are associated with supernovae — resulted from the intensive study of 
a relatively small number of bursts without regard for whether these bursts formed a well 
defined statistical sample. However, realizing that for any one burst we are only sampling 
an anisotropic radiation pattern from one direction, we now have to study distributions of 
burst properties to reconstruct what bursts look like from all directions. Well-defined burst 
samples are required to derive these distributions. 
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Unfortunately, we do not yet have a large sample of bursts with redshifts, which in 
most cases are required for calculating the intrinsic burst properties. The two dozen or 
so bursts with redshifts were detected by various detectors with different detection thresh- 
olds and energy sensitivities, and the follow-ups that resulted in the redshifts depended 
on the vagaries of weather and telescope availability. BATSE provided a large burst sam- 
ple with well-understood thresholds, but without direct redshift determinations for most of 
these bursts. However, redshifts can be determined indirectly from the lag-luminosity or 
variability-luminosity relations. Here we calculate self-consistently the redshifts for a large 
subset of the BATSE bursts using the lag-luminosity relation. To accommodate the prob- 
lematic burst GRB980425, which is significantly underluminous if the lag-luminosity relation 
is a single power law, Salmonson (2001) and Norris (2002) proposed that the relation should 
be a broken power law. This reduces the luminosity of many long-lag bursts, moving this 
population closer to the observer. 

We Use the resulting redshifts to calculate both the peak bolometric luminosity Lb and 
the isotropic energy E iso ; both quantities assume that the burst radiates isotropically. We 
consider the isotropic energy E iso to be more fundamental, and therefore model the energy 
distribution. 

The energy distribution can discriminate between the two proposed models for the 
structure of the jets. The uniform jet model (Frail et al. 2001; Bloom, Frail, & Kulkarni 
2003) assumes that all jets have a constant surface energy density e (energy per solid angle) 
but differing opening angles 6 0 . The model makes no predictions for the energy distribution. 
In the universal jet profile model (Rossi, Lazzatti, & Rees 2002; Zhang & Meszaros 2002) 
all jets have the same surface energy density e as a function of off-axis angle 9, and thus 
the observed differences in Lb or E 1S0 results from our angle 9 V relative to the jet axis. This 
model predicts that if the surface energy density e is a power law with index k, e oc 9 k , then 
the energy distribution will be a power law with index a = 1 — 2/k, where N(Ei SO ) oc E~% ; 
hence a = 2 for k = — 2. A Gaussian profile results in a — 1. Numerical modelling of 
hypernovae (Zhang, Woosley, & MacFadyen 2003) show that the surface energy density of 
the outflows should vary with off-axis angle, but do not predict the profile. Thus the profiles 
are empirical models to fit the data. 

The simple universal jet profile model can be falsified by the observed energy distribution 
but the uniform jet model cannot. However, physically plausible refinements, such as allowing 
the parameters of the jet profile to vary somewhat from burst to burst — the jet profile is now 
only quasi-universal — results in much greater flexibility in reconciling predicted and observed 
energy distributions (Lloyd- Ronning et al. 2003). Nonetheless, we compare our calculated 
energy distributions to the predictions. 
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In these calculations we use a cosmology with H 0 = 70 km s -1 Mpc -1 , f l m = 0.3 and 
Qa = 0.7. We use the notation p(a | b) to mean the probability of a given b. 

In §2 we present the methodology used in this study: calculating the burst redshifts 
(§2.1); fitting the intensity distribution (§2.2); and using the cumulative probability to test 
the quality of the fit (§2.3). The implementation of this methodology is discussed in §3 and 
the results are provided by §4. Finally, our conclusions are in §5. 


2. Methodology 
2.1. Calculating Redshifts 

In the absence of a large number of measured spectroscopic redshifts, we use the lag- 
luminosity relation (Norris et al. 2000) to calculate peak luminosities Lb, and subsequently 
redshifts, for a large burst sample. Assume that for a set of bursts we have measured the lag 
between the lightcurve in two energy bands r 0 (s), the peak photon flux P (ph s _1 cm -2 ) in 
the energy band E 1 -E 2 , and the spectral fit for the time period over which P is measured. 
The spectral fit is characterized by the parameters a , /? and E p of the “Band” model (Band 
et al. 1993). The instrument team is assumed to have converted the observed peak count 
rate into the peak photon flux using the fitted spectral model and the response function. 

The lag tb in the burst’s frame and the observed lag To, both between the same energy 
bands in their respective frames, are related by time dilation (which increases the lag) and 
spectral redshifting (which shifts the smaller lag at high energy into the observed band). 
The first effect is simple and universal — a factor of (1 + z) -1 . The second depends on the 
burst’s spectral evolution, and may vary from burst to burst. Thus 

tb = T o9(z) • ( 1 ) 

We will assume that g(z ) is universal, and as an ansatz we assume g(z) = (1 + z) C t . Time 
dilation contributes — 1 to cv, while the redshifting of temporal structure with a smaller lag 
from higher energy contributes a positive constant (e.g., ~ 1/3; Fenimore et al. 1995 found 
~ 0.4). 

Empirically the lag has been related to the apparent bolometric peak luminosity Lb 
(erg s _1 ; Norris et al. 2000) 

L b = Q(tb ) • (2) 

The observed peak bolometric energy flux Fb is related to the bolometric peak lumi- 
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nosity L b through the redshift, and an assumed cosmology: 

F b = L b /[47tD 2 l ] , (3) 

where D B is the luminosity distance. Note that L B is the “isotropic” peak luminosity, the 
peak luminosity if the observed flux were beamed in all directions. If the flux is actually 
beamed into a solid angle Afl, then the actual peak luminosity is only Af2/4-7r of L B . We 
define 


(E) = 

f b /p 

( 4 ) 




where P = 

/ N(E)dE 
J Bi 

( 5 ) 

and F b — 

poo 

/ EN(E)dE ; 

Jo 

(6) 


N(E) is the photon spectrum. 

The implicit equation for redshift, 

p Q(r 0 g(z)) 

4mD\(E) ’ { } 

must be solved for each burst. The functional form of the lag- luminosity relation Q can be 
calibrated from the small set of bursts for which z is known (see below). We have an ansatz 
for g(z) with one unknown constant c T . This constant can be calculated using the dependence 
of t on energy bands (Fenimore et al. 1995). Alternatively c T can be fit simultaneously with 
the determination of Q. 

For this large database we then have L B , z, To, P and the spectral parameters for the 
peak of the lightcurve. This database can be used to study the distribution of intensity 
measures such as the total burst energy and the peak luminosity. 


2.2. The Intensity Distribution 

We use a Bayesian analysis to investigate intensity distributions. In this study we 
consider the total isotropic energy E lso as the relevant intensity measure; the isotropic energy 
is the energy emitted if the observed flux were radiated isotropically and not beamed towards 
us. Note that the same formulae apply substituting L B for E iS0 . 

Assume our burst database provides the isotropic energy E- lso , the peak flux P and the 
redshift z for each burst (Band 2001). The threshold peak flux P m ; n is assumed to be known, 
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from which the threshold isotropic energy £'i SO ,min for each burst can be calculated: to first 
order -Bi SO ,min = E lS o Pmm/P- The energy distribution function is the normalized probability 
distribution p{E\ so \ a v Mj,z, I) where aj- is the set of parameters that characterize the jth 
model distribution function represented by Mj, and I specifies general assumptions about 
the distribution function. For generality the energy distribution function has been given a 
redshift dependence. However, the observed E lso are not drawn from p(Ls | a}, Mj, z, I) but 


from 


P{.E\so ■^'iso,min i Mj , Z, /) 


p(E lso cij , Mj , z, I)9(E lso E \ so 

,min) 

ijS is0 niin < ^^ ,so P(E'lSO 1 Mj, Z, I ) 


( 8 ) 


where 9{x) is the Heaviside function (1 for positive x and 0 for negative x); E^ is drawn 
from the observable part of the luminosity function. We now form the likelihood 


N 

A? = _|_ j[_ p{P'iso,i | £'iso,min,ii Mj, Zj, /) 
i=l 


(9) 


where the zth burst has isotropic energy E^, threshold isotropic energy E iso , m i n ,i and redshift 

z i * 

In the “frequentist” framework best-fit parameters are typically found by maximizing 
A j. The Bayesian analysis is based on p(dj\D,Mj,I), the posterior probability for the 
parameters; here D is the set of observed peak bolometric luminosities and their thresholds. 
Note that A j — p(D | d}, Mj, I ). Thus an estimate of the parameters is 


(aj) 


J daj aj p[aj \D,Mj, I ) 


/ daj aj p(D | aj, M j} I)p(aj \ M jt I) 
f daj p(D | aj,Mj , /)p(dj \M 5 ,I) 


J* daj Oj p(E ls0 j 

-£^so,min,i ? aj,Mj,z i ,I)p(a j \Mj,J) 

f daj [jj_j p(£'iso,i ^iso.minA? 


( 10 ) 


where p(dj \ Mj, I ) is the prior for the parameters a}. If E j = p(D | a}, Mj, I)p(a 3 \ Mj, I ) is 
sharply peaked, then the expectation value of the parameters occurs at the peak of Ej. Note 
that Ej is the likelihood in eq. 9 times the priors for the parameters, and is proportional to 
the posterior probability p(dj \D,Mj,I). The posterior probability is also used to determine 
the acceptable parameter range, which is often more meaningful than the “best” parameters 
values. 
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2.3. The Cumulative Probability 

The Bayesian approach does not provide a goodness-of-fit statistic. However a frequen- 
tist statistic can be derived. For each burst the cumulative probability is 

poo 

P(E'\so,\ I -E'iso ,min,i ? Qj 5 i M j} I)= p(E iso \ ■^iso,min,i5 5 Zi, Mj, I) dE iso . (11) 

J E\so,i 

If the assumed energy distribution function is an acceptable characterization of the obser- 
vations (which would be the case if the model Mj is correct) and all the assumptions are 
valid (e.g., the cosmological model is correct), then the cumulative probabilities P(E lS0ti ) 
should be uniformly distributed between 0 and 1, and have an average value of (P(E is0ti )) = 
1/2 ± (12 iV) -1 / 2 for N bursts in the sample. 


3. Implementation 
3.1. Datasets 

Our data are all derived from various BATSE burst databases. We start with a database 
for 1438 bursts that includes the lags and their uncertainties, the peak flux P on the 256 ms 
timescale, the burst duration T 90 , and hardness ratios among the 4 BATSE energy channels 
(30-50, 50-100, 100-300, and 300-2000 keV). Of these 1438 bursts, 1218 have positive lags. 

To calculate the average energy (E) (see eq. 4) we use the parameters of “Band” spec- 
trum fits by Mallozzi et al. (1998) to the 16 channel BATSE “CONT” count spectra accu- 
mulated over the peak flux time interval (usually 2.048 s) for the 580 bursts in our database 
that are also in the Mallozzi et al. database. For the 858 bursts that are not in the Mallozzi 
et al. database we assume the average low and high energy spectral indices of a = —0.8 
and /? = —2.3 found by Preece et al. (2000). Because E p , the energy of the peak of the 
EN(E) oc vf„ curve, and HR32, the 100-300 keV to 50-100 keV hardness ratio, are strongly 
correlated, we use the empirical relation E v = 240 HR32 keV (for HR 32 < 2.25) for the bursts 
without spectral fits. Figure 1 shows the resulting scatter plot for E p vs. HR32. 


3.2. The Lag-Luminosity Relation 

Based on less than a dozen bursts, Norris et al. (2000) found that L B oc r^ 1 ' 15 . However, 
this single component relationship predicts a much greater L B than observed for GRB980425, 
the burst that appears to coincide with the supernova SN1998bw. By introducing a break 
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in the lag-luminosity relation Salmonson (2001) and Norris (2002) were able to include 
GRB980425. The resulting lag-luminosity relation is: 

L- ol = 2.18(r/0.35 s) Ci , c L = -1.15 for r < 0.35 s, - 4.7 for r > 0.35 s. (12) 

where L 51 = Lb/ 10 51 erg s -1 . In our calculations we use both the single (i.e., cl = —1.15 
for all r) and two component lag-luminosity relations. 


4. Results 

4.1. Lag-Luminosity Relation and Redshift 

We first investigate the lag-luminosity relation. Salmonson (2001) and Norris (2002) 
introduced a break in the simple power law relation to include the long lag, very low lumi- 
nosity GRB980425. As a result of this break other long lag bursts in the BATSE database are 
assigned low luminosities, and consequently small distances; Norris suggests that there may 
be a population of nearby low luminosity bursts. Figure 2 demonstrates this: the solid lin e 
is the cumulative distribution of burst redshifts assuming a single power law lag-luminosity 
relation, while the dashed line shows the distribution for the broken power law. Introducing 
a break in the power law shifts the low redshift bursts closer. Note this break is included only 
to include GRB980425 in the lag-luminosity relation; if GRB980425 was not associated with 
SN1998bw or if GRB980425 was anomalous, then this break is unnecessary and unsupported 
by other data. 

Figure 2 also shows that there are few bursts with z > 10. These bursts may indeed 
be at high redshift; finding such high redshift bursts is one of the goals of the Swift mission. 
However, on theoretical grounds bursts with z > 17 are not expected. Given the dispersion 
in the lag-luminosity relation, it is not surprising that some bursts are assigned redshifts as 
high as z = 65; large redshifts result from high luminosities for bursts with very small lags, 
which are particularly difficult to measure accurately. Thus the lag-luminosity relation does 
not give unphysical results at the high luminosity end. Norris (2002) placed an upper limit 
on the luminosity when he found a large number of very high redshift bursts. 

We use the simple single power law lag-luminosity relation without any cutoffs or limits 
in the following analysis. 

For our database the median redshift is z m = 1.58. As can be seen from Figure 1, 
this redshift is greater than the redshift range where the form of the lag-luminosity relation 
makes a difference. 
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4.2. The E lso Distribution 

We focus on the distribution functions of the bolometric energy. To perform quantitative 
estimates of distributions as a function of redshift (e.g., of the burst rate) the threshold for 
the database must be known, and therefore studies have often focused on the peak luminosity 
Lb] L b is more closely related to the peak flux P, which has a relatively sharp instrumental 
threshold P m in . However, to study distributions of intrinsic burst quantities the threshold 
for that quantity for each observed burst will suffice. Thus our database is sufficient to study 
both the isotropic energy E iSO and the peak luminosity Lb distributions. 

Figure 4 shows the scatter plot of the isotropic energy E iso . The region just above the 
threshold is underpopulated, suggesting that the thresholds are underestimated. Therefore, 
we raise the threshold to P m ; n = 0.5 ph cm -2 s _1 . This raises the median redshift from 1.58 
to 1.62. 


4-2.1. Power Law E 1S0 Distribution 

Universal jet structure models predict a power law energy distribution function with an 
index of a = 2 (p(E iso ) oc E T“); for a finite number of bursts the distribution must have 
a low energy cutoff E 2 . We therefore fit our burst sample with this distribution function. 
Figure 5 shows the resulting likelihood surface which peaks at a = 1.76 and log(P 2 ) = 50.2. 
This value of E 2 is the smallest value of E iso in our sample; burst samples that extend to 
fainter bursts will probably have smaller observed values of E lso , and consequently E 2 is most 
likely smaller. The contour plot indicates that the determination of oc in independent of the 
value of E 2 . The 95width of 0.05. 

The likelihood alone does not indicate whether the power law model is a good descrip- 
tion of energy distribution. The average of the cumulative probability (§2.3) should be 
( P(E iS0 )) = 1/2 ± (12N)~ 1 / 2 for N bursts. Figure 6 shows that the cumulative P(E iSO ) for 
P m in = 0.5 ph cm -2 s -1 is very close to the straight line that is expected if the distribution 
function is a good description of the data. Although ( P(E 1S0 )) = 0.4642 ±0.0089 (N = 1054) 
formally differs by 4 a from the expected value of 1/2, considering the likely systematic uncer- 
tainties (e.g., in the value of P m j n , in the determination of the lags, in the corrections for the 
redshifting of high energy lightcurves), this value of ( P(E lso )) indicates that the distribution 
can be represented by a power law. 
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4-2.2. Lognormal Distribution 

The second model distribution we consider is a lognormal energy distribution (see §2.5.1 
of Band 2001); the parameters are the centroid E lsocen and the logarithmic width a E of the 
distribution. For this distribution the posterior distribution (relevant to a Bayesian analysis) 
as a function of In E^ cea and a E and the likelihood (relevant to a frequentist analysis) are 
the same if the prior is constant in In Fiso.cen and o E . Figure 7 shows the likelihood surface; 
the peak occurs at E isooeD = 3 x 10 51 ergs and o E = 2.7. However, as the figure shows, a 
broader distribution (larger a E ) with a smaller centroid E lso>cen is not ruled out. 

For Pmin = 0.5 ph cm -2 s -1 (P(E tso )) = 0.4821 ±0.0089 ( N = 1054), which is consistent 
with (P(Eiso)) = 1/2 at the 2a level. Note that for P m j n = 0.3 ph cm -2 s -1 (P(E iso )) = 
0.4598 ± 0.0085 ( N = 1162), which is differs from the expected value by nearly 5a. Figure 8 
shows the distribution of the cumulative probability for P m i n = 0.5 ph cm -2 s _1 ; as can seen, 
the observed cumulative probability distribution is consistent with the expected distribution. 


5. Summary 

We use the lag-luminosity relation to calculate self-consistently the redshifts for 1218 
BATSE bursts. For the bursts without the spectral parameters required by the calculation 
we use average low and high energy spectral indices, and a peak energy E p derived from the 
hardness ratio. We find that the redshift is quite sensitive to the spectral parameters. We 
use both single power law and broken power law lag-luminosity relations, and find that the 
broken power law relation does indeed predict a population of low luminosity, nearby bursts. 
For both forms of the relation the median redshift is 1.58. 

We use the redshifts to calculate the apparent peak bolometric flux and the isotropic 
energy, both assuming that the bursts radiate isotropically. 

We find that our burst data can be fit by a power law energy distribution with a = 
1.76 ± 0.05 (95% confidence); considering only the statistical uncertainty the power law dis- 
tribution is formally not a good fit, but with the likely systematic uncertainty the power law 
distribution is probably a good description of the data. While our power law fit is inconsis- 
tent with the original universal jet profile model (which predicts a = 2), it is consistent with 
the jet profile models where parameters are permitted to vary. 

A log-normal energy distribution also describes the data; the data permit a smaller 
average energy if the distribution is wider. More faint bursts will bound the lower end of 
this distribution. 
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Redshift 

Fig. 2. — Cumulative distribution of bursts by redshift. Solid line — one component lag- 
luminosity relation. Dashed line — two component relation. Both distributions have the 
same median redshift of 1.6. 
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Fig. 3. — Bolometric peak luminosity Lb vs. redshift z for single component lag-luminosity 
relation. The lines represent the threshold peak luminosity for a threshold peak flux of 
P > 0.3 ph cm -2 s -1 and E p — 30 (2 dots-dashed line), 100 (dot-dashed), 300 (solid) and 
1000 keV (dashed). 
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Fig. 4. — Scatter plot of E lso vs. the i?; SO| min- On the dotted line E iso = E iso>min . The 
paucity of bursts just above the dotted line suggests that P min should be raised from 0.3 to 
0.5 ph cm -2 s -1 . 
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Fig. 5. Contour plot of the likelihood function for a power law energy distribution. A 
threshold peak flux of P min = 0.5 ph cm -2 s' 1 was imposed. The power law has an index of 
a and a low energy cutoff of E 2 . 
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Fig. 7. — Contour plot of the likelihood function for a lognormal energy distribution. A 
threshold peak flux of P min = 0.5 ph cm" 2 s" 1 was imposed. The distribution is characterized 
by a centroid energy and a logarithmic width. 






